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Abstract. Consider the Navier-Stokes flow past a rotating obstacle with a general time- 
dependent angular velocity and a time-dependent outflow condition at inflnity. After rewrit- 
ing the problem on a flxed domain, one obtains a non-autonomous system of equations with 
unbounded drift terms. It is shown that the solution to a model problem in the whole space 
case R'' is governed by a strongly continuous evolution system on L%(R'^) for 1 < p < oo. 
The strategy is to derive a representation formula, similar to the one known in the case of 
non-autonomous Ornstein-Uhlenbeck equations. This explicit formula allows to prove U'-L'' 
estimates and gradient estimates for the evolution system. These results are key ingredients 
to obtain (local) mild solutions to the full nonlinear problem by a version of Kato's iteration 
scheme. 

1. Introduction and main result 

In this paper we consider a model problem in M"^ for the flow of an incompressible, viscous 
fluid past a rotating obstacle with an additional time-dependent outflow condition at infinity. 
The equations describing this problem are the Navier-Stokes equations in an exterior domain 
varying in time with an additional condition for the velocity field at infinity. 

In order to motivate our model problem, let O C M*^ be a compact obstacle with smooth 
boundary, let Q := M'^ \ O be the exterior of the obstacle and let m G C([0, oo); M''^'^) be a 
continuous matrix-valued function. Then, the exterior of the rotated obstacle at time t > is 
represented by Vt{t) := Q{t)Q where Q{t) solves the ordinary differential equation 

r dtQ{t) = m{t)Q{t), t>0, 

[ Q{0) = Id. ^^'^^ 

With a prescribed velocity field Voo G ^^^([0, oo);]^'^) at infinity, the equations for the fluid on 
the time-dependent domain r2(t) with no-slip boundary condition take the form 

vt- Av + v ■Vv + Vq = in fi(t) x (0, oo), 

divv = in il.{t) x (0, oo), 

v{t, y) = m{t)y on d9.{t) x (0, oo), (1.2) 

lim v{t, y) = Voo{t) for t € (0, oo), 

ly|->oo 

v{0,y)=uoiy) in fi, 

where v and q are the unknown velocity field and the pressure of the fluid, respectively. 
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The disadvantage of this description is the variabihty of the domain 0,{t), and the fact that 
the equations do not fit into the L^-setting, due the velocity condition at infinity. Assume 
for the time beeing that m{t) is skew symmetric for t > 0; this implies that for all t > the 
matrix Q{t) is orthogonal. Then, by setting 

x = Q{tfy, u{t,x) = Q{tf{v{t,y)-Voo{t)), p{t, x) = q{t,y), (1.3) 

the above equations can be transformed to the reference domain Q and the new velocity field 
u vanishes at infinity. Then ()1.2p is equivalent to the following system of equations 

Ut- Au- M{t)x ■\7u + M{t)u ^ 

+Q{t)'^Voo{t) ■ Vu - Q{tfdtVoo{t) > = in X (0, oo), 

+?x • Vu + Vp J 

divii = in X (0, oo), 

u{t,x) = M{t)x - Q{t)'^Voo{t) on 517 X (0,oo), (1.4) 
lim u{t,x) = for t G (0,oo), 

|a;|— >oo 

m(0, x) = uo{x) in $7, 

where M{t) := Q{t)^m{t)Q{t). The main difficulty in dealing with this problem arises since 
the term Ai{t)x ■ V has unbounded coefficients. In particular, the lower order terms cannot 
be treated by classical perturbation theory for the Stokes operator. 

Note that even if we assume that m{t) = m is independent of time (this implies that 
also A4{t) = Ai is independent of time), equation (|1.4p is still non-autonomous due to the 
time-dependent first order term Qit)"^ Voo ■ V (except in some special cases discussed below). 

However, by using localization techniques similar to |GHH06| . this problem is finally reduced 
to a model problem in M"^ and a model problem in a bounded domain. Since Q{t)dtVoo{t) = 
F{t), t > 0, i.e. it is constant in space, we may put this term in the pressure p. Hence, in this 
paper we discuss the following linearized model problem in M'^ 

ut-Au- {M{t)x + fit)) ■ Vu + M{t)u Vp = in M*^ x (0, oo), 

divu = in M'^ X (0,oo), (1.5) 

■u(O) = uq in M'^, 

where we allow gener al coefficients M € C([0, oo); M'^^^) and / G C([0, oo); M*^). If we set 
M{t) := Q{t)'^m{t)Q{t) and f{t) := —Q{t)'^Voo{t) then we obtain the linearization of equation 
()1.4p with = M'^. Such a model problem also arises in the analysis of a rotating body with 
translational velocity —Voo{t), see |Far05| . 

Existence and uniqueness of a mild solution of an autonomous variant of problem (jl.2p 
without an outflow condition, i.e. Voo = 0, and m{t) = m, was investigated in quite a few 
papers, see |His99a| . |His99b| . |GHH06| and |HS05| . Hishida was even able to deal with a time 
dependent rotation in |His01| . however only for angular velocities of a special form. 

For the problem including an additional outflow condition at infinity, there are only a 
few results. Indeed, in the special case, where m{t)x = io{t) x x and u : [0, oo) is 
the angular velocity of the obstacle and v^o '■ [0, oo) ^ a time-dependent outfiow velocity, 
Borchers |Bor92| constructed weak non-stationary solutions for the equations (|1.4p . Moreover, 
Shibata |Shi08| studied the special case where m{t) = m, foo(*) = 'Voo and mvoo = 0. The 
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condition mvoo = 0, i.e. Q{t)'^Voo = kvoo for k G { — 1, 1}, ensures that ()1.4p is still an 
autonomous equation and the solution of (|1.4p is governed by a Co-semigroup which is not 
analytic. The physical meaning of the additional condition mvoo = is that the outflow 
direction of the fluid is parallel to the axis of rotation of the obstacle. The stationary problem 
of this latter situation was analysed in |Far05| . 

The assumption mVoo = was recentely relaxed by the second author in |HanlO| . Indeed, 
he was able to deal with the model problem in M'^ where m{t)voo 7^ and Voo{t) = Voo- 
However he assumes that m{t) and m{s) commute for all t,s > which can physically be 
interpreted by the fact that the axis of rotation is fixed. 

The aim of this work is to remove the latter additional condition, i.e. m{t) and m{s) need 
not to commute and Voo may be time-dependent. 

As usual the Helmholtz projection P allows us to rewrite (|1.5p as an abstract Cauchy 
problem in La{W^), where L§(M'^) denotes the space of all solenoidal vector fields in LP(R)^: 



Note that it immediately follows from |HS05| that for fixed t > 0, the operator A{t) is the 
generater of a Co-semigroup, which is not analytic. The fact that the semigroup is not analytic 
prevents us from employing standard generation results for evolution systems, see |Paz83| 
Chapter 5] and references therein. For the same reason, L^-L'^ estimates and gradient estimates 
don't follow from standard arguments. 

Therefore, we first derive a representation formula for the solution of (jl.Sp . In order to 
derive this representation formula we transform (jl.Sp to a non-autonomous heat equation 
which can be explicitly solved, see Section [3j It turns out that the transformation to a non- 
autonomous heat equation is crucial to deal with our problem in this generality since the 
different transformation used in |HanlO| caused the additional assumption that M{t) and 
M{s) commute for all t,s > 0. 

In the following we denote by {f/(i, s)}t,s>o the evolution system on generated by the 
family of matrices {—M{t)}t>o, i-c 



Note that dsU{t,s) = U{t,s)M{s). 

We are now ready to present our main result. 

Theorem 1.1. Letl < p < oo, M € C([0, oo); M'^'^'^) and f € C([0, oo);M'^). The the solution 
of (|1.6p is governed by a strongly continuous evolution system {T(t, s)}t>s>o C C{L^{M.'^)'^). 
Moreover, the evolution system {T(t, s)}t>s>o admits the following properties: 



u'{t) - A{t)u{t) =0, t > 0, 
u{Q) = uq. 



(1.6) 



Here: 



A{t)u 
D{A{t)) 



P (An + {M{t)x + f{t)) ■ Vn + M{t)u) 

{u G W^'P{m.y n Ll{W^) : M{t)x ■ Vu G LP{Ry}. 




(1.7) 
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(a) For To > set A/tq := sup{||f/(i, s) || : t,s G [0,To]}. Then for 1 < p < oo and 
p < q < oo there exists C := C{Mto, d) > such that for u € L^{W^) 



1 1 



\\T{t,s)u\y^(^d^<C{t-s) Hp W||n||iP(ig.), 0<s<t<To, (1.8) 

||Vr(t,s)n||^,(K.) <C(t-s)-^(^~^)~^||n||iP(Kd), 0<s<t<ro. (1.9) 

/n particular, if the evolution system {U{t,s)}s,t>o is uniformly bounded, i.e. Mtq < 
M , for some M > and all Tq > 0, we may set Tq = oo. 
(b) For I < p < q < oo, s > and u € L^iW^) we have 



d 1 1 
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lim (t-s)nP -^J\\T{t,s)u\\Li^^d^=0 and lim {t - s)2\\VT{t,s)u\\LP^^d) = 0. 

t^s, t>s ^ ' t^s, t>s ^ ' 

Next we consider the nonlinear problem 

n'(t)-A(t)n(i)+P((n(t)-V)t/(t)) = 0, t>0, ^^^^^ 

n(0) = tio, 

with initial value uq G Lct(M'^). 

For given < Tq < oo, we cah a function u G C([0, Tq); Lg(]R'^)) a mz/d solution of (fLTn]l if 
satisfies the integral equation 

u{t) =T{t,0)uo - [ T(t,s)P((n(s) • V)n(s))ds, t > 0, (1.11) 
Jo 

in Lct(M'^). By adjusting Kato's iteration scheme (see |Kat84) ) to our situation the existence 
of a unique (local) mild solution follows, cf. |HanlO| for details. 

Corollary 1.2. Let 2 < d < p < q < oo, M e C([0, oo); M'^^'^), / G C([0, oo); M*^) and 
uq G La{M.'^). Then there exists Tq > and a unique mild solution u G C([0, Tq); Lct(M'^)) of 
il.lO\) . which has the properties 

t-^{-v~du{t) G C([0,ro);L^(]R"')), (1.12) 

if(|-i)+lvu(t) G C{[0,To);L'^{Ry'^). (1.13) 

If p < q, then in addition 

t^(^"^)||u(t)||i,(Kd)+t^||Vn(t)||iP(Kd) ^0 ast^O. (1.14) 

Moreover, in the case d = p we may set Tq = +oo provided \\uo\\i^d(^d^ is small enough and 
{U{t, s)}s^t>o is uniformly bounded. 

Remark 1.3. In particular, {Uit., s)}s,t>o is uniformly bounded if M{t) is skew symmetric 
for alH > 0. 



2. Proof of Theorem 11.11 

Let M be as in Theorem 11.11 and let {U{t,s)}s,t>o be the evolution system on that 
satisfies ()1.7p . We consider the system of parabolic equations of the form 

{ dtu{t,x) - A{t)u{t,x) = 0, t>s,xeW^, 

\ u{s,x) = (p{x), X G M'^, 
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for s > fixed, initial value ip G LP{W^)'^ and some p G (l,oo). Here the family of operators 
A{t) is of the form 

A{t)u{x) := (^Aui{t, x) + (M(t)x + f{t), Vui{t, x))^ ^ - M{t)u{t, x), t>0, x eM'^. 

As in |GL081 Lemma 3.2] or |HanlO| . we first develop an explicit representation formula. To 
be more precise, we show in Section [3] that for p G (1, oo) and ip G L^(]R'^)'^ the solution u to 
()2.ip is governed by a strongly continuous evolution system {T(t, s)}t>s C C{D-'{W^)'^) which 
is explicitly given by 

u{t, x) : = (r(t, s)ip){x) := {k{t, s, •) * ip){U{s, t)x + g{t, s)), t > s, x gW^, (2.2) 



where 



kit, s, X) := (4,)./2(dLQ,.)i/^ ^^^' .)e-i<«M^.^)dy, t > . > 0, x G (2.3) 
g{t, s) := f U{s, r)f{r)dr, Qt,s := f U{s, r)U*{s, r)dr, t > s > 0. 

J S J S 

Similar to |DPL07| one can show that for ip G C~(M'^)'^ the solution u of ([23]) given by ([22]) 
is a classical solution. 

A simple calculation shows that divT'(t, s)(p = for ip G C^(M'^) and t > s > 0. Hence, the 
restriction T{t,s) := T(t, s)\]^P(^^d^ is an evolution system on La{W^)- In particular, u{t) : = 
T(t, 0)no is a solution to (|1.6p . 

By similar arguments as in the proofs of |GL08| Lemma 3.2] or |Hanl01 Lemma 2.4], for 
To > there exists C := C{d,MTQ) > (see Theorem 11.11 for the definition of Mj-g) such that 



1 



\\Q,J\\<C{t-s)-^, 0<s<t<n, ^^^^ 

(detQt,s)^ > C(t- s)i, 0<s<t<ro. 

Moreover, if is uniformly bounded in Tq we may write Tq = oo in (j2.4p . 

Proof of Theorem We start by showing the estimate (jl.Sp . Let Tq > 0. By the change of 
variables ^ = U{s, t)x and by Young's inequality we obtain 

\\T{t,s)u\\Li(^^d) < \detU{s,t)\i\\k{t,s,-)\\Lr^^d)\\u\\LP^^d), t> s>0, 

where 1 < r < cx) with ^ + ^ = 1 + |- Further, by the change of variable y = Q]/gZ we obtain 



l-r 



mt,s,-WLr^^d) = \mt,s)\\ [ -i-,e"^ (detQt,,)^dz 

JMrf \(47r)2 J 

< C7||C/(t,s)||(detQt,s)^, t>s>0, 



for some C > 0. Now ([M]) yields (fLS]) . 

To prove the gradient estimate (|1.9p . we first observe that 



VT(t, s)u(x) = j u{U{s, t)x + g{t, s)k{t, s, y) {U'^{s, t)Q~ly)'^ dy, t> s>0, xe M'^. 
Now, ()1.9p follows similarly as above. 
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Since (|2.ip is uniquely solvable for ip € C'^iW^)'^^ see Section [3l the law of evolution is valid, 
i.e. 

f{t,s)^ = f{t,r)f{r,s)^, (2.5) 

holds for < s < r < t and every 99 G C~(M'^)°'. The density of C^{W^Y in U'{W^Y yields 

that ([23]) even holds for all tp G L'P{W^f. 

In order to prove the strong continuity of the map (t, s) 1— )• r(t, s) on < s < t we apply 

1/2 

the change of the variables y = Q^\. z, to see that 

f{t, sMx) = -^U{t, s) ■ [ ip{U{s, t)x + g{t, s) - Q|,z)e-^dz 
(47r) 2 JRd 

holds. For t > s fixed, we pick two sequences (tn)nGN and {sn)n&n such that t„ > s„ holds for 
every n G N and (t„, s^) — )• (t, s) as n — >■ 00. For every ip G C^(M'^)'^ and every x G M'^ we 
now obtain 

tf{U {Sn, tn)x + g{tn, S„) - Qt„,s„z) (p{U (s, t)x + C/(t, s) - Ql^z) 

as n — )■ 00. Lebegue's theorem now yields T{tn, Sn)^ — >■ T(t,s)ip as n — >■ 00 for every (p G 
C^(]R'')'^. The density of C^{W^y in LP(M'')'' implies the strong continuity. 

In order to prove Theorem [LU b) let u G L^{R'^), t-s < 1 and choose (Mn)nGN C C^(]R'^) C 
La{M.'^), such that lim„_!.oo ll'W~^n||LP(R'*) — 0- The triangle inequality together with the L^-L'^ 
estimates (|1.8p imply that there exist constants Ci , C2 > such that 

(t_s)^(^-i)||r(t,s)n||i.(jj,) 

< (t _ 5)1 (f-i) ||T(t, S){U - Un)h%m + - S)^^^"') \\T{t, s)Un\\Ll 



<Ci\\u-Un\\LP,^d)+C2it- s)^yP "^WUnhll'^d), < t - S < 1, n G N. 



Hence, lim(t — s)2Vp iJ \\T{t, s)u\\i^i md\ = by letting first t ^ s and then n — ?> oo. The 



second assertion in Theorem 11.1( b) is proved in a similar way. □ 

3. Representation Formula 



In this section the representation formula ()2.2p is derived. The general idea is to do a 
coordinate transformation in order to eliminate the unbounded drift and the zero order term 
of the operator A{t). For this purpose we set 

z := U {s, t)x + g{t, s), 

where 

g{t,s) := / U{s,r)f{r)dr, 

J s 

and we look for a solution u of (j2.ip with initial value ip G C^{M.'^)'^ in the form 



u{t, x) = U {t, s)w{t, U{s, t)x + g{t, s)). 



(3.1) 
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By recalling (|1.7p we obtain from a straightforward computation that 

dtu{t, x) = -M{t)U{t, s)w{t, z) + U{t, s)(^{U{s, t)M{t)x + C/(s, Vwi{t, z))) ^ 

+ C/(t,s)9j'u;(t,z), 
holds. Moreover, we can write equation p.ip component-wise as 

Ui{t, x) = C/jj (t, s)wj{t, U{s, t)x + s)), for i = 1, . . . , d, 
i=i 

and thus for the spatial derivatives of u we obtain 

d 

Vui{t,x) =^Uij(t,s)U*{s,t)Vwj{t,z), 
i=i 
d 

V^Ui{t,x) = ^Uij{t,s)U*{s,t)V^Wj{t,z)U{s,t). 
i=i 

In particular, the drift term can be written as 

d 

{M{t)x + fit), Vuiit, x)) = U^j{t, s){U{s, t)M{t)x + U{s, t)f{t), Vwj{t, z)). 

i=i 

Thus, the function u solves problem (|2.ip if and only if for every i = 1, . . . ,d, the function 
Wi : — >■ M is a solution to 

( dtWi{t,z) = TT[U{s,t)U*{s,t)V'^Wi{t,z)], t>s,zeM.'^, 
\ Wi{s,z) = ipi{z), z € M"'. 

By our transformation we now obtained an uncoupled system of parabolic equations with 
coefficients only depending on t. More precisely, for i = l,...,d, the equation ()3.2p is a 
non-autonomous heat equation. It is well known that such a problem can be uniquely solved 
(cf. |DPL071 Proposition 2.1]) and that for every ipi G C^(M'^) its unique solution is explicitly 
given by the formula 

w,it,z) = ^ [ ^,{z - y)e-^<«*"'^'J'>dy, (3.3) 

(47r)2(detQt,,)2 Jr^ 

where 

Qt,s = J^U{s,r)U*{s,r)dr. (3.4) 

Now, via (|3.ip . the unique solution to our original problem (j2.ip is given by the representation 
formula 

u{t, x) = {k{t, s, •) * u){U{s, t)x + g{t, s)), (3.5) 

where the kernel k{t,s,x) is defined in (|2.3p . 

Note that the right hand side of (|3.5p is even well defined for each L^(]R'^)'^-function ip. Thus, 
this explicit formula can be used to define an evolution system on LP{R'^f in the following 
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way. For G WiW^)'^ we set 

if for t = s, 

s, x) * ^){U (s, + (7(t, s)) for t > s. 



Since problem (|3.2p is uniquely solvable it follows via (|3.ip that T{t, s)ip is the unique 



oo fmd\d 



solution of (j2.ip for initial value ip £ 
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